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Temperature dependence of the electron spin g factor in CdTe and InP 
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Temperature dependence of the electron spin g factors in bulk CdTe and InP are calculated and 
compared with experiment. It is assumed that the only modification of the band structure related 
to temperature is a dilatation change in the fundamental energy gap. The dilatation changes of 
fundamental gaps are calculated for both materials using available experimental data. Computations 
of the band structures in the presence of a magnetic field are carried out employing five- level Pp 
model appropriate for medium-gap semiconductors. In particular, the model takes into account spin 
splitting due to bulk inversion asymmetry (BIA) of the materials. The resulting theoretical effective 
masses and g factors increase with electron energy due to band nonparabolicity. Average g values 
are calculated summing over populated Landau and spin levels properly accounting for the thermal 
distribution of electrons in the band. It is shown that the spin splitting due to BIA in the presence 
of magnetic field gives observable contributions to g values. Our calculations are in good agreement 
with experiment in the temperature range of K to 300 K for CdTe and K to 180 K for InP. The 
temperature dependence of g is stronger in CdTe than in InP due to different signs of the band- 
edge g values in the two materials. Good agreement between the theory and experiment strongly 
indicates that the temperature dependence of spin g factors is correctly explained. In addition, we 
discuss formulas for the energy dependence of spin g factor due to band nonparabolicity, which are 
liable to misinterpretation. 

PACS numbers: 72.25.Fe, 71.70.Ej, 72.25.Rb, 78.47.-p 



I. INTRODUCTION 

Temperature dependence of the bulk electron spin g 
factor in semiconductors is of interest both for scientific 
reasons as well as for possible spintronic applications. As 
far as the theory is concerned, a correct description of 
electron spin properties can test validity of the k-p the- 
ory for nonzero temperatures. As to the experiment, it 
is now possible to measure the electron spin g value up 
to the room temperature using quantum beats in a time- 
resolved photoluminescence of spin states and related ef- 
fects. The latter furnished a consistent experimental pic- 
ture of g(T) for electrons in GaAs (see Ref. 1 and the ref- 
erences therein). It was shown that one can successfully 
describe the g factor if, in the nonparabolic k-p the- 
ory one takes the dilatation change of the energy gap, 
in agreement with previous theoretical predictions 2-4 . 
On the other hand, one should account for the obvious 
fact that, as the temperature increases, more and more 
Landau and spin levels are populated by electrons. In 
nonparabolic conduction bands of III-V compounds the 
g factors change with electron energy from their values 
at the band edge to the free electron value +2 at high 
energies 5,6 . Since g depend on electron energy, one mea- 
sures in reality their values averaged over all populated 
levels. 

In our present paper we are concerned with tempera- 
ture dependences of the bulk electron g factors in CdTe 
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and InP. The temperature variations of g in these ma- 
terials have been measured but the theoretical descrip- 
tions are missing 7,8 . Ito et al 8 attributed almost all the 
temperature variation of g in CdTe to the far-band con- 
tributions which was clearly not justified, as the authors 
themselves recognized. Thus, the subject can be regarded 
as controversial, as was previously the case for GaAs, see 
Refs. 7, 9, 10. 

The band structure of CdTe and InP is similar to that 
of GaAs. In CdTe the g value is negative at the band 
edge and it tends to +2 going through zero as the energy 
increases. In InP the g value is positive at the band edge 
so that, tending to +2, its energy dependence is much 
slower. Thus we deal with two distinctly different cases. 
Another aspect of the present work is the spin splitting 
due to the bulk inversion asymmetry (BIA), called the 
Dresselhaus splitting 11 . This splitting leads to a number 
of effects, also in the presence of an external magnetic 
field. In our description we account for the spin splitting- 
due to BIA at finite magnetic fields and discuss this con- 
tribution. To our knowledge, this problem has not been 
considered before. 

Our work has three objectives. First, we describe the 
temperature dependence of the electron g factors in CdTe 
and InP and compare the theory with existing experi- 
mental data. Second, we analyze the temperature depen- 
dences of g values in the two cases and show that they 
result from the opposite signs of g at the band edges 
of both materials. Third, we study the effect of Dres- 
selhaus spin splitting at finite magnetic fields on the g 
value. Finally, we discuss validity of a frequently em- 
ployed formula for the energy dependence of g factor in 
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III-V and II- VI compounds and indicate how it should 
be used. Our paper is organized in the following way. In 
Section II we summerize the band structure calculations 
and indicate how the average g values are computed. In 
Section III we describe our calculations for CdTe and 
compare them with experimental data, Section IV con- 
tains similar program for InP. In Section V we discuss 
our results and in Section VI we summerize them. In 
Appendix we consider the energy dependence of g value, 
liable to misinterpretation. 



II. 



THEORY 



InP and CdTe are medium-gap semiconductors (MGS) 
and a three-level k-p description, successfully used for 
narrow gap semiconductors 12,13 , is not adequate for de- 
scribing their band structures. The reason is that in MGS 
the fundamental gap E between the Tg and Tg levels is 
not much smaller than the gap E\ between the Y% level 
and the upper r 7 conduction level. It has been demon- 
strated that an adequate way to treat the conduction 
band of MGS is to use a five- level model (5LM), which 
is equivalent to 14 bands (including spin) in the k-p de- 
scription (see Refs. 14, 15 and the references therein). 
According to the five-level model the spin g value at the 
conduction band edge is [14], 
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where E Po = 2m a P 2 /h 2 , E Pl = 2m P 2 /h 2 , G = E +A a 
and Gi = E\ + Ai. The spin-orbit energies Ao and Ai 
relate to (r^, T|) and (Ty, r§) levels, respectively, A is 
the interband matrix element of the spin orbit interac- 
tion between the (T^, rg) and (r 7 , r§) multiplets (see 
[14, 16]), and C is due to far-band contributions. Fig- 
ure 1 shows schematically the five-level k-p model used 
in our calculations. For A = Eq. (1) reduces to the 
formula given first by Hermann and Weisbuch 17 . Cal- 
culating the electron energies away from the band edge 
one deals with the effects of band's nonparabolicity and 
inversion asymmetry. In particular, an appearance of the 
matrix element Q is due to the bulk inversion asymme- 
try. Since the 5LM for electrons in the presence of a 
magnetic field and its use for magnetooptical properties 
of MGS was described in some details before 14 ' 15 , we only 
mention here the main elements of this approach. Thus 
the model includes exactly the Ty, Tg, Tg, Vj and r§ lev- 
els at the center of the Brillouin zone and the resulting 
k-p matrix has dimensions 14 x 14. There exist three 
nonvanishing interband matrix elements of momentum: 
Po, Pi and Q. If one takes Q = and k z — (where hk z 
is the momentum along the magnetic field) the 14 x 14 
initial matrix factorizes into two 7x7 matrices for the 



spin-up and spin-down states. These matrices are soluble 
by envelope functions in the form of harmonic oscillator 
functions and the eigenenergy problem for different Lan- 
dau levels (LLs) n reduces to diagonalization of 7 x 7 
determinants. 

Taking the asymmetric gauge for the vector potential 
A = [-By, 0, 0] one obtains the harmonic oscillator func- 
tions in the form exp(ik x x + ik z z)<& n [(y — yo)/L], where 
ya = k x L 2 . Here L = (h/eB) 1 / 2 is the magnetic radius. 
If the Q element is included (it comes from an inversion 
asymmetry of MGS crystals) the initial 14 x 14 matrix 
does not factorize and is not soluble in terms of a single 
column of harmonic oscillator functions. Physically, this 
means that the resulting energy bands are not spherical. 
Since the nonsphericity of the conduction bands in MGS 
is small, one can find the eigenenergies looking for the en- 
velope functions in terms of sums of harmonic oscillator 
functions (see Ref. 18). This leads to number determi- 
nants composed of the fundamental 7x7 blocks on the 
diagonal coupled by nondiagonal parts involving Q and 
k z elements. The eigenenergies are computed truncating 
the resulting big determinants. In our computations we 
used typically 112 x 112 determinants. All calculations 
were performed taking a magnetic field B parallel to the 
[001] direction. 

Next, we consider average values of the spin g factor 
measured as a function of temperature. The measure- 
ments are usually done in relatively pure samples hav- 
ing low free electron densities. The electrons are excited 
across the gap into the conduction band and into both 
spin states. The spin states are almost equally popu- 
lated and the circularly polarized light produces a well 
defined coherence between them. The excited electrons 
quickly thermalize and are distributed among Landau 
levels (LLs) according to the lattice temperature with- 
out losing their spin or phase. Then they interfere and 
quantum beats in the photoluminescence or other effects 
are observed from many LLs. According to this picture 
the observed signal represents an average over the popu- 
lated levels. The electron thermal distribution over LLs 
determines their contribution to the average g value. 

We assume the fc z -dependence of electron energies in a 
simplified form 
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where n is the LL number, ± signs correspond to the 
two spin states, k z is the wavevector along the direction 
of B, and TOq is the effective mass at the band edge. 
The description of energies contains the intricacies of 
the band structure mentioned above. The spin g value is 
defined as (in formulas we use g* symbol) 



g* = (£+-£- k )/f, B B. 



(3) 



An averaging procedure involves a summation over n 
and integrations over k x and k z . A simple calculation 
gives the average value of g* in the form 
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FIG. 1: Five-level Pp model of band structure in medium- 
gap semiconductors CdTe and InP. The zero of energy is cho- 
sen at the F§ edge. Interband matrix elements of momentum 
Po, Pi, Q as well as interband matrix element of spin-orbit 
interaction A are indicated, C and C' symbolize far-band con- 
tributions to the band-edge effective mass and spin g factor, 
respectively. 



¥{T) = ^ , (4) 

where 

A= S4^^ rd£ ' (5) 

and 

in which the summation is over the LLs, f(£ , £) is the 
Fermi-Dirac distribution function, and the square roots 
come from the integrations over k z . The integrations 
begin from the lower of the two states E l n for each n, 
which can be either £+ or £~ depending on the sign of 
the g value. 

The average g value, as given by Eq. (4), is affected 
by the temperature in two opposite ways. As the tem- 
perature T increases and the absolute value of the funda- 
mental gap Eq decreases, the spin g$ value at the band 
edge decreases. On the other hand, with increasing T the 
electrons populate higher LLs and band's nonparabolic- 
ity comes more and more into play. The latter is known 
to make the g value less negative (see Refs. 5, 6). Thus, 
as T increases, the average g* decreases or increases de- 
pending on the relative strength of the two effects. We 
emphasize that we do not use in our calculations Eq. (1), 
it is quoted only to make clear the dependence of </q on 
Eq and other parameters. 



III. CdTe 
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FIG. 2: Dilatation gap Eq(T) in CdTe versus temperature 
calculated from Eq.(7) with the use of experimental values of 
D, dEo/dP and a th (T), see Eq. (2). 



As mentioned above, it was demonstrated that the 
temperature change in the effective mass and the spin 
g factor in a material is governed by a dilatational vari- 
ation of energy gaps, of which the fundamental gap is of 
primary importance. Thus one needs to determine the 
dilatational change in the fundamental gap due to tem- 
perature since the directly measured total temperature 
change is due to both the dilatation of the crystal lattice 
and its vibrations (phonons). We quote the determina- 
tion of Eq 1 (T) for CdTe since, to our knowledge, it has 
not been carried out before. The dilatational change in 
the gap is given by 

AE il (T) = -3D(^j^j\ th (T')dT' , (7) 

where D is the bulk modulus, 8E /dP is the pressure- 
induced gap shift, and a t / l (T) is the linear thermal ex- 
pansion coefficient (see also Ref. 19). The quantities D 
and 8E /dP are directly measurable, for CdTe there is 
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FIG. 3: Calculated spin g factor in bulk CdTe versus electron 
energy in the conduction band. Full points are g values for 
consecutive LLs and k z = 0, solid line indicates g value for 
n = as a function of k z , empty points and dashed line: the 
same as above but neglecting bulk inversion asymmetry ( i.e. 
putting Q = 0). 



dE /dP = 0.08 eV/GPa 20 ' 21 (see also Refs. 22-26) and 
D = 42 GPa 26 . The linear thermal expansion coefficient 
a t h{T) was calculated 28 and measured 29 . Using these re- 
sults we performed the integration indicated in Eq. (7) 
and obtained the dilatational variation of the fundamen- 
tal gap in CdTe shown in Fig. 2. It is seen that E|f (T) 
goes through a maximum at T « 65 K but, all in all, the 
temperature variation between K and 300 K is rather 
small. 

Knowing Ef(T) and assuming that other gaps, mo- 
mentum matrix elements and far-bands contributions do 
not depend on the temperature, one can perform the 
band structure calculations outlined above if the band 
parameters are known. For CdTe we take the following 
parameter values at T = 0: E Pa = 21.07 eV, E Pl = 5.1 
eV, Eq = 13.29 eV (see Ref. 30), E = -1.6 cV 31 , A = 
-0.95 eV 32 < 33 , E 1 = 3.76 eV 34 < 35 , Ai = 0.27 eV 30 , A = 
-0.19 eV 30 , C = -0.5355 and C'=-0.0129. The far-band 
contributions are taken to obtain at T = the band edge 
values mj = 0.093 mo 36 and g$ = -1.66 9 . The zero of en- 
ergy is chosen at the Tg edge, see Fig. 1, so the energies 
above are positive while the energies below are negative. 

Figure 3 shows the result of our intermediate calcula- 
tions for CdTe, given as example. The calculations are 
performed for fixed values of T and B. In order to in- 
vestigate the effect of BIA on the g value at finite B we 
carried out the computations in two versions: 1) using the 
full 5LM, i. e. including the matrix element Q which, as 
mentioned above, incorporates BIA (full points and solid 
line); 2) putting Q = 0, i. e. neglecting the effect of 
BIA (empty points and dashed line). The full points in 
Fig. 3 indicate calculated spin g values for consecutive 
LLs (at k z = 0) beginning with LL n = 0. It is seen that 
the g value increases with the LL number n (or, equiv- 



alently, the energy) due to band's nonparabolicity. The 
solid line shows the calculated g for n = as a function 
of k z (or, equivalently, the energy). It can be seen that, 
for lower LL numbers n, the g value behaves very simi- 
larly for n-dependence and k z - dependence. As argued by 
the present authors 37 , the above dependences should be 
identical within the description by the three-level model, 
see also Appendix. The empty points and dashed line 
show the corresponding quantities calculated with Q = 
0. Here both g*(n,k z — 0) and g*(n — 0,k z ) are prac- 
tically the same. The differences between the full points 
(and solid line) and the empty points (and dashed line) 
are directly due to the effect of BIA. 

It is well known that, in the absence of magnetic field, 
the spin splitting due to BIA is AS = 2^\k 2 {k 2 x k 2 + 
fc 2 fc 2 + kyk z ) - Qklkykl] 1 / 2 . For the presence of B we 
can not give an analytical expression for the spin split- 
ting due to BIA. However, one can qualitatively say that, 
at B ^ 0, k x and k y components are replaced by (nB) 1 ' 2 
or [(n + V^B] 1 / 2 terms, while the k z component along the 
magnetic field remains the same. For n = the "trans- 
verse" components are small (or zero) and for S « the 
longitudinal component fc 2 w 0, so that BIA gives almost 
no contribution, which agrees with the results shown in 
Fig. 3. The increasing energy S corresponds to the in- 
crease of n or k z (or both), so the contribution of BIA 
grows. This is reflected by an increasing difference of the 
results for Q ^ and Q = shown in Fig. 3. 

We follow the similarity of g*(n =const, k z ) and g*(n, 
k z = const) dependences shown in Fig. 3 assuming g*(n 
=const, k z ) to be equal to to g*(n, k z — const) in the 
summation of LLs and integration over k z , see Eqs. (5) 
and (6). This approximation is in fact quite good since 
the region of high energies: S > 90 meV for CdTe (see 
Fig. 3), where this approximation begins to break down, 
is weakly populated by electrons and it gives only small 
contribution to the average g value. 

Using the above assumption we performed the summa- 
tion and integration indicated in Eqs. (5) and (6) and 
determined the average value of g according to Eq. (4). 
Our final results for CdTe are shown in Fig. 4. It can 
be seen that the experimental values of Oestreich et al 7 
(lower temperatures) and that of Ito et al 8 (room temper- 
ature) are described very well. For temperatures above 
200 K the theory is slightly higher then experiment. We 
conclude that the main reason behind the observed in- 
crease of the spin g value with temperature is that, as 
the temperature increases, higher Landau and spin levels 
contribute to the cverage g value. Since the band struc- 
ture predicts an increase of g with growing energy (see 
also the discussion in Appendix) the resulting average g 
becomes higher (less negative). This mechanism is well 
illustrated in the inset of Fig. 4. 
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FIG. 4: Spin g factor in bulk CdTe versus temperature. 
Empty points: experimental data of Ref. 7, full point: data 
of Ref. 8. Solid line: theoretical average g factor calculated 
according to Eq. (4). Inset shows how consecutive LLs con- 
tribute to the average g value at T = 300 K. 



InP is a medium-gap semiconductor similar in many 
respects to GaAs. However, the electron spin g factor 
at the band edge of InP is positive, in contrast to InSb, 
InAs, GaSb and GaAs 17 . Since, as we mentioned above, 
at high electron energies g tends to +2 (if one neglects 
the effects of BIA) there remains not much room for the 
energy variation of g value. 

The dilatational change in the fundamental energy gap 
of InP was estimated by Hazama et al 4 . However, we 
revise this estimation because a part of the procedure 
adopted in Ref. 4 was based on calculations, while in our 
approach we use exclusively experimental information. 
To determine the dilatation gap of InP we use again for- 
mula (7). We take D = 71 GPa 38 < 39 , and dE /dP =0.084 
eV/GPa 40 . As to the function a t h{T), it was measured 
by various authors, see Refs. 41-43. We follow the mea- 
surements of Haruna et al 41 , which basically agree with 
those of other authors but are more complete. In the 
data of Ref. 41 we correct the point at T = 8 K (since 
the given value ceth(8K) has the wrong sign) by interpo- 
lating between the measured values at K and 15 K. The 
complete function ath{T) is then used for the numerical 
integration indicated in Eq. (7). Our final results for 
Eq 1 (T) are shown in Fig. 5; the obtained variation in 
the dilatation gap is noticeably smaller than that given 
in Ref. 4. We use the results indicated in Fig. 5 in our 
further procedure. 

To perform the band structure calculations we take 
the following band parameters 44 : Ep = 20.93 eV, Ep 1 
= 0.165 eV, E Q = 15.56 eV, E = -1.423 eV, A = -0.108 
eV, Et = 3.297 eV, A x = 0.201 eV, A = 0.08733 eV 45 , C 
= -2.467 and C = -0.08045. The far-band contributions 
are taken to obtain at T = the band edge values 
= 0.07927m 36 and g^ = +1.204. Our band structure 
calculations are performed similarly to those described 
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FIG. 5: The same as in Fig. 2 but for InP. 



above for CdTe. 

In Fig. 6 we show intermediate calculated results for 
g*(n,k z = 0) and g*{n = 0,k z ) according to the 5LM 
in two versions: 1) including the matrix element Q, i.e. 
including BIA (full points and solid line); 2) putting Q 
= 0, i.e. neglecting BIA (empty points and dashed line). 
Again, as discussed above, in the region of energies £ < 
190 meV the two dependences are very similar (for Q 
— they practically coincide for all energies) so that, 
when performing the integration over k z in Eq. (5) we 
assume that the fc z -dependence of g value for each LL is 
the same as that given by the black points. As far as the 
contribution of BIA to the spin splitting is concerned, 
the picture is similar to that for CdTe: at low energies 
the contribution of BIA vanishes but it grows with the 
energy. Clearly, in our final calculations we do take into 
account the effect of BIA. 

Figure 7 shows our final results for the temperature 
dependence of the average g value in InP, computed ac- 
cording to Eqs.(4)-(6). It is seen that the available exper- 
imental data of Oestreich et al [7] are described very well. 
Similarly to CdTe, the main reason for the increase of g 
value with the temperature is that, as T grows, more 
Landau and spin levels in the nonparabolic conduction 
band of InP are populated with electrons. On the other 
hand, the temperature variation of average g in InP is 
considerably weaker than that in CdTe for reasons given 
above. 



V. DISCUSSION 

One should bear in mind that the g values measured by 
various authors do not agree too well with each other. For 
CdTe, Meyer et al 46 used the spin resonance to measure 
g values in the temperature range of 4.2 K < T < 66 
K and obtained the data which can be described by the 
linear dependence g*(T) = -1.682 + 2.97xl0~ 4 T. All 
these values are somewhat lower than those shown in 
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FIG. 6: The same as in Fig. 3 but for InP. 
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FIG. 7: Spin g factor in bulk InP versus temperature. Empty 
points: experimental data of Ref. 7. Solid line: theoretical 
average g factor calculated according to Eq. (4). Inset shows 
how consecutive LLs contribute to the average g value at T = 
300 K. 



Fig. 4. Sprinzl et al 47 used the spin quantum beats to 
measure the g values that agree with those shown in Fig. 
4 at low temperatures but arc considerably lower near 
room temperature. We believe that the data shown in 
Fig. 4 are reliable since the lower temperature results of 
Oestreich et al 7 are consistent with the result of Ito et 
al 8 for T = 300 K. 

As for InP, we note that the low-temperature result of 
Weisbuch and Hermann 48 g'£ — +1.26, obtained with the 
use of spin resonance, is distinctly higher than the val- 
ues shown in Fig. 7. However, such differences between 
measured g values are not unusual, see Fig. 4 of Ref. I 
for GaAs. We were unable to find experimental measure- 
ments of the g factor in InP for temperatures above 200 
K. 

As follows from Figs. 3 and 6, the g values calculated 
taking Q ^ 0, i.e. using the complete 5LM description, 



are for higher electron energies considerably higher than 
those obtained taking Q = 0. This means that the agree- 
ments between experiment and theory shown in Figs. 4 
and 7 give clear indications that the influence of bulk in- 
version asymmetry on the spin splitting at B + 1 is not 
negligible in CdTe and InP. This evidence is still some- 
what indirect since the average values are given by sums 
over many LLs and integrations over k z . However, as 
each term of the sum over LLs and the integration over 
k z is affected by BIA, the evidence for its non- negligible 
effect seems convincing. One could verify more directly 
the effect of BIA in the presence of a magnetic field by 
investigating at low temperatures an anisotropy of the 
spin splitting with respect to the orientation of a mag- 
netic field. It should be mentioned that the effect of BIA 
on the g values of heavy holes in GaAs quantum wells 
was recently investigated by Kubisa et al 49 . 

Finally, we believe that the distinctly different tem- 
perature dependences of spin g values in CdTe and InP, 
both experimental and theoretical, are significant. The 
increase of g in CdTe between T « and T = 1 80 K is 
Ag* « 0.16, while in InP the increase in the same tem- 
perature range is Ag* « 0.049. As we remarked above, 
this difference is due to different energy dependences of g 
factors in the two materials because at high energies the 
g factor should reach the free-electron value of +2 (if one 
neglects BIA), while the initial band-edge value in CdTe 
is #0 « -1.66 and in InP it is g$ w +1.20. In consequence, 
g* {£) dependence in InP is distinctly weaker than that in 
CdTe. The above reasoning is confirmed experimentally 
which strongly indicates that our interpretation, as de- 
scribed above, is correct. We emphasize this conclusion 
since the reason for the temperature dependence of the 
spin g factor remains to be a matter of controversy, cf. 
Ref. 10. 



VI. SUMMARY 

Dilatational changes in the fundamental energy gaps 
of CdTe and InP are determined from available experi- 
mental data in order to use them in calculations of band 
structures at nonzero temperatures. The five-level Pp 
model is employed to compute the band structures of 
these medium-gap semiconductors in the presence of a 
magnetic field. In particular, energy dependences of the 
electronic spin g factors due to band nonparabolicities are 
obtained. Next, average g values are calculated for differ- 
ent temperatures summing over populated Landau and 
spin levels and integrating over longitudinal momentum 
hk z . The increase of <?'s with growing temperature are al- 
most exclusively due to the population of higher Landau 
and spin levels in nonparabolic conduction bands. The 
calculated g*(T) dependences are compared with avail- 
able experimental values and good agreement between 
experiment and theory is obtained for both materials. 
The temperature increase of spin g factor is stronger in 
CdTe than in InP, which is related to the negative band- 
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edge value of g in CdTe and the positive one in InP. It is 
shown that the bulk inversion asymmetry gives observ- 
able contributions to the spin splittings in the presence 
of a magnetic field. Frequently used formulas for the 
energy-dependent spin g factors in nonparabolic conduc- 
tion bands of III-V compounds are discussed. 



Appendix 

Since our work is closely related to dependence of the 
electron spin g factor on the electron energy £ in III- 
V and some II- VI semiconductor compounds, we discuss 
here the use of an often employed formula for g*{£). To 
the best of our knowledge, the first description of g*(£ ) in 
narrow-gap III-V compounds was given in 1961 by Lax, 
Mavroides, Zeiger, and Keyes 5 (LMZK) and a different 
but equivalent formula was derived in 1963 by Zawadzki 6 . 
The formula of LMZK, although old, is still used in the 
literature (see for example Refs. 8, 50-52) so it merits a 
discussion. 

It is important to see how the formula is derived. The 
underlying formulation is due to Bowers and Yafet 12 , see 
also Ref. 6. The three-level P p description (see Refs. 
12, 53) takes into account the Tg conduction level and 
Tg, valence levels. The resulting P-p 8x8 Hamilto- 
nian is solved in terms of harmonic oscillator functions 
neglecting small free-electron terms. Final equation for 
the energies is obtained in the form 



£(£ - E )(£ - E - A )+ 



The additive terms +1 in Eq. (A.3) and +2 in Eq. (A.4) 
result from the free-electron contributions in the P-p 
theory, which were neglected in the original Bowers and 
Yafet treatment 12 , cf. Ref. 14. Equations (A. 2), (A.3), 
(A.4) amount to the LMZK formulas in which band's 
nonparabolicity enters via the energy dependence of m* 
and g* (LMZK formulas do not contain the free electron 
terms). 

The problem with the above formulation is that the 
energy £^ k on the LHS of Eq. (A. 2) is the same as 
that entering m*(£^ k ) and g*(£^ k ) on the RHS. Sup- 
pose one fixes B, n and k z and tries to calculate the spin 
splitting. It is clear that it is not enough to take the en- 
ergy difference for the plus and minus signs in Eq. (A. 2) 
and calculate the energy on LHS because both the mass 
m*(£^ k ) and the spin factor g*{£^k ) have different val- 
ues for different signs of the spin. In particular, if n or 
k z are large, the orbital term in Eq. (A. 2) is consid- 
erably different for the opposite spins which should be 
necessarily taken into account. In fact, the form given 
by LMZK and used in the literature is inconvenient since 
the energy £^ k should be determined self-consistently. 
It is considerably simpler to use directly Eq. (A.l) and 
solve for the energies. However, in most applications one 
simply employs formula (A.4) for the energy dependence 
of g* (£) forgetting the orbital term and the fact that in 
a nonparabolic band g*(£^k ) i s different for each spin 
direction. 

Finally, we note that a formulation analogous to that 
given by LMZK for 3LM has been derived by the present 
authors for GaAs-type medium gap semiconductors with 
the use of 5LM, see Refs. 37, 54. The result is that 



-PM^+l)+kl](£-E - 2 -A )T-P^A oS = , (A.l) Eq ' is stiU vaM but the formulas for m *^J and 



9*( £ nk z ) become 



where s = eB/h and other symbols have been de- 
fined above. For specified B,n,k z and the spin sign ±, 
Eq. (A.l) represents a cubic equation for the energy 
£(n, fcz,±). In order to proceed further, we divide Eq. 
(A.l) by {£ — Eq){£ — E — A ) and after a simple alge- 
braic manipulation obtain for the conduction band 



to 



1+C- 



+ E_ 



Pi 



1 

+ w 



+ 



AA^E^E^ / 1 



HeB 



^nk z 



where 



1, 



21.2 



K Z k 



to 



E 



Pa 



2 1 



in which 



= 2 + 



EZ 



G 



Eq = E - £„ kz 



G = E a + A - £ nk ^ 



(A.2) 
(A.3) 

(A.4) 

(A.5) 
(A.6) 



and 



E t G E G t 



1 



4Ay/Ep E Pl 



&0 <JQ 



1 



E 1 G E U ± 



(A.7) 

1 1 

Gf-Et 

(A.8) 



where Eq and Gq are still given by Eqs. (A.5) and (A.6), 
respectively, and 



E t = E l~£n 



Gf = E 1 + A 1 -£± kz 



(A.9) 



(A.10) 
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The terms proportional to Ep 1 are related to the inter- 
band Tg-r| and F^-Fy interactions, see Fig. 1. On both 
sides of Eq. (A. 2) the energy is E nk and, for given B, 
n, k z and ±, this energy should be determined in a self- 



consistent way. The above formulas do not contain the 
element Q, i.e. they neglect the effects of bulk inversion 
asymmetry. 
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